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Abstract 

We investigate the production and possible detection of gravitational waves stemming from 
the electroweak phase transition in the early universe in models of minimal walking technicolor. 
In particular we discuss the two possible scenarios in which one has only one electroweak phase 
transition and the case in which the technicolor dynamics allows for multiple phase transitions. 
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I. INTRODUCTION 



Recent progress in the understanding of the phase diagram of generic asymptotically 
free gauge theories (l]-[6l has led to a renewed interest in this class of models GHH. F° r a 
recent review of the latest developments see [1J. Explicit examples of technicolor models, 
not in conflict with electro weak precision tests, have been put forward in ||2Tl4ll9l[T0|. The 
simplest incarnations of these models are known as (Ultra) Minimal Walking Technicolor 
models l|2l44l [10l[TT| and indicated in short by MWT and UMT respectively The principal 
feature is that the gauge dynamics is such that one achieves (near) conformal dynamics 
for a small number of flavors and colors. It has been shown that one can construct 
cold dark matter candidates via either the lightest technibaryon, here termed Technicolor 
Interacting Massive Particles (TIMP)s |10l - [T8| , or new heavy leptons naturally associated 
to the technicolor theory ||T91 - |2T| . The TIMP is naturally of asymmetric dark matter 
type JUrO/ meaning that its relic density does not have a thermal origin. Within the 
(U)MWT models such a relic density has been estimated in I1T01 ITTIIT5I . Weak isotriplet 
TIMPs (iTIMP)s have been shown to be interesting candidates of dark matter in [|22[|. 

Another interesting cosmological arena is the temperature driven electroweak phase 
transition within technicolor theories |23l - l2"6l| . We have investigated in much detail this 
phase transition for the MWT and UMT models using the low effective Lagrangian ap- 
proach. We discovered that there is a sizable region of the low energy effective theories' 
parameters yielding a sufficiently strong first order electroweak phase transition to drive, 
in principle, electroweak baryogenesis. We have also discovered quite a rich phase dia- 
gram in the case of the UMT model ||2"5H and more generally whenever several underlying 
matter representations are simultaneously present in the technicolor dynamics [24J. An 
interesting problem is if such a transition is observable via detecting the cosmological 
gravitational waves (GW)s produced at the transition itself. We will describe the topic 
of GWs in more detail in the next section. Whether or not these waves are observable 
depends on the strength of the electroweak phase transition. We will investigate this 
issue using the two concrete models discussed above. To be able to study the production 
of GWs, we need to use a slightly improved treatment of the phase transition compared 
to our earlier work. We confirm the results of ||2~3l |25l , and find that the MWT model 
Lagrangian can support a sufficiently strong electroweak phase transition leading to an 
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observable signal at BBO [|27|. An interesting feature of the UMT model is the presence of 
multiple electroweak phase transitions arising at different temperatures. This is due to the 
interplay between two distinct chiral phase transitions, one directly responsible for the 
electroweak symmetry breaking and the other decoupled from the standard model (SM). 
However we find that the planned experiments searching for GWs will have hard times 
discovering the signal originating from this model and more sensitive ones are needed. 

It is, however, possible to increase the strength of the first order phase transition by 
considering partially gauged technicolor models ||3l H|. They have, by construction, a 
large number of techniflavors but only two of them are gauged under the electroweak 
symmetry. This choice reduces the contribution to the electroweak precision parameters, 
while the large number of techniflavors enhances the strength of the first order phase 
transition. The nonrenormalizable axial anomalous contributions to the effective low 
energy potential is partially responsible for increasing the strength of the transition when 
increasing the number of techniflavors. We will, however, investigate the spectrum of 
gravitational waves associated to the electroweak phase transition stemming from generic 
models of partially gauged technicolor elsewhere. For other simple models which predict 
potentially strong GWs, see, for example, ||28ll29| . 

Summarizing we investigate in detail the MWT and UMT models, at the effective 
Lagrangian level, and show that the MWT can lead to detectable gravitational waves 
while the UMT cannot. 

II. GRAVITATIONAL WAVES PRODUCTION SETUP 

In this section we lay the basics of the GW production from strong first order phase tran- 
sitions and we present the relevant parameters of our theories needed for the calculation 
of the gravitational signal. 

First of all, let us review how a first order phase transition takes place in the early 
universe and why it can produce GWs. In a first order phase transition there are two 
distinct minima separated by a potential barrier. The phase transition can be thought 
to start taking place at the moment where the two vacua are at the same energy level. 
Immediately after, the true-vacuum state to be lowers its potential level compared to 
the one of the false vacuum and therefore despite the existence of the barrier, quantum 
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mechanically, there is a finite probability for the system to pass from the false to the 
true vacuum. The phase transition occurs through nucleation of bubbles of the true 
vacuum ||30l . (For a nice exposition see [31].) The nucleation is possible due to quantum 
tunneling and thermal fluctuations. The bubbles of the true vacuum expand until they 
cover the whole space, which means that the phase transition has been concluded. For the 
production of GWs a quadrupole moment is required, and since the bubbles are spherical, 
and therefore have no quadrupole, it seems at first sight that no production can take place. 
However, there are at least two different ways of producing GWs from bubbles. The first 
one is when bubbles collide. Apparently in this case the spherical symmetry is destroyed 
and GWs are produced. The second source of GW production is due to turbulence of the 
plasma because of the bubble's motion. 

For a given theory, there are basically two parameters that determine the signal of the 
GWs produced Il32l - l34l . The first one a is defined as the ratio of the latent heat e of the 
phase transition at the bubble nucleation temperature over the energy density of the false 
vacuum. Practically the latent heat is the energy released as the system tunnels from the 
false to the true vacuum. It is given by 

dV 

e = -AV - TAs = -AV + T— , (1) 

oT 

where V is the potential. Since we are interested in theories where the phase transition 
takes place around the electroweak scale (~ 250 GeV), the energy density is dominated by 
the radiation part. The parameter a practically measures how strong the phase transition 
is. As we shall discuss later on, large a, i.e., strong first order phase transition, leads 
to enhanced amplitude for the GW and therefore better detectability. From this point of 
view, theories with strong phase transitions are more interesting. 

If roughly speaking a affects the amplitude of the GWs, the second model dependent 
parameter jS, determines the characteristic frequency. This is because corresponds to 
the rate of change of the nucleation probability and therefore has units of inverse time. 
This means that is approximately the duration of the phase transition, and provided 
we know the velocity of the bubble expansion, it determines the size of the bubble (having 
ignored the initial size which is negligible compared to the final). Let us see this explicitly. 
The bubble nucleation rate at nonzero temperature is given by 



where Se = S3/T and 



S 3 = 



drAnr 




+ Vty, T) 



(3) 



is the Euclidean 3-dimensional action, (p is the bubble profile. We are looking for a 
least action solution which has an 0(3) symmetry. The equation of motion exhibiting 
manifestly this symmetry and dictated by the minimization of the Euclidean action is 



where r is the radial coordinate. The boundary conditions are d(p(r = 0)/dr = 0, i.e., we 
require the solution to be smooth at the center of the created bubble, and <p(r = 00) = 0, 
meaning that far away from the bubble the system is still in the false vacuum. If one 
imagines <p to be the position of a particle and r to be the time, the above equation 
corresponds to the equation of motion of a particle within a potential - V((p) with unit 
mass and a Stokes type of drag force proportional to the velocity given by the second 
term of the equation above. In general, the bounce solution cannot be found analytically 
due to the complexity of the equation. However within the thin wall approximation a 
closed solution has been found. The thin wall approximation is valid when the difference 
in the height of the two minima is small compared to the barrier. Let us call (p e the "escape 
point", i.e., the value of <p(r = 0) of the solution of Eq. Q. In the particular case where 
the two minima are almost degenerate, <p e should be very close to the true minimum to 
reduce the work done against the drag force: The particle starts from a nearly flat point in 
the potential and therefore it would take some time (r in this particular case) to build up 
its velocity, and then go downhill fast in order to come at rest again at the false minimum. 
Practically, this means that within this approximation, the friction term of Eq. Q can be 
safely ignored and an analytical result can be obtained. It also means that the change 
from one minimum to the other happens "fast" and the profile of <p is quite sharp. This 
justifies the name "thin wall." Generally, if one is not sure whether or not the thin wall 
approximation is valid, a numerical solution of Eq. Q is needed. The standard way of 
finding this numerical solution (which is the one we also used in order to get our results) 
is by guessing the value of the <p e . We know that (p e should be between the two minima. 
If our guessed value of (p e is closer to the true minimum that in reality, solving Eq. Q 



with the boundary conditions (p(r = 0) = <pg and d(p(r = 0)/dr = will "overshoot 
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[ (p 2d(p dV((p,T) 



= 0, 



(4) 



dr 2 



r dr d(p 



the solution, meaning that after some r the solution will start going to -co. On the other 
hand, if (/)' nal is closer to the false minimum than the actual <p e , the trial solution will 
"undershoot" the real solution meaning that it will never reach the false vacuum. 

Having introduced the Euclidean action, /3 is defined as the time derivative of the 
action —dS^/dt, evaluated at the nucleation temperature that we shall introduce shortly. 
In the early universe the expansion parameter a ~ T" 1 and therefore the Hubble parameter 
H = (l/a)da/dt = -(l/T)dT/dt. Consequently 

j6 ^dS E ^d(S 3 /T) 

It is understood that everything is evaluated at the nucleation temperature. This tem- 
perature is defined as the temperature where the rate of bubble nucleation per Hubble 
volume and time is approximately one. This means 

r^H 4 -» Tin— * -— , (6) 

ntpi 4 

where m P i is the Planck mass and we used H - T 2 /m P i. The above equation gives the 
bubble nucleation temperature. We shall denote this temperature by T* and the value of 
the Hubble parameter at nucleation by H» below. Recall that S3 is known once we have 
found the bounce solution, and substitute it in Eq. ([3]). 

The parameters a and (} are the essential input parameters we need from the specific 
model under investigation. The strength and the frequency of the gravitational signal 
produced by the first order phase transition are encoded in these two parameters. Let us 
review the basic arguments of how GWs are produced due to bubble collisions in a more 
quantitative way, following the scaling argument presented in [|32|. GWs are produced 
through quadrupole (or higher moment) emission. For the quadrupole, the GW power is 
P = (G/5){(Qj.) 2 ), where G is the Newton constant and Q i; is the quadrupole moment of 
the transverse and traceless part of the energy-momentum tensor. Note the dependence 
of the power on the triple derivative of the quadrupole with respect to time, something 
which is also true in electromagnetism. The quadrupole moment has dimensions of mass 
times distance squared and therefore dimension analysis dictates that the triple derivative 
would have the units of kinetic energy over time. Not all the energy gained from tunneling 
from the false to the true vacuum is in the form of kinetic energy. If k is the fraction of 
the latent heat in the form of kinetic energy (the rest being heat), E kin ~ ta*PradOWj8) 3 - 



As we have already mentioned the latent heat is ap ra d, and we also have multiplied by 
the volume of the bubble ~ {v^lff where V\, is the velocity of the bubble walls. From 
Friedmann's equation we know that H 2 ~ Gp CT n, and therefore we can trade G for H*. In 
addition, p crit = (1 + «)p ra d- Using all of the above and keeping in mind that E G w = P/fi, 
we get that Qgw ~ (H*lfi) 2 k 2 a 2 v 3 b l(l + a) 2 . This is as good as dimension analysis can get 
us. 

The GW production due to bubble collisions was first studied in f|33l |35] - |37| . These 
calculations were based on numerical simulations of bubbles colliding using the so-called 
envelope approximation, which consists of considering only the nonoverlapping regions 
of the collided bubbles as sources of GW production. In such case [33J 

\ jS I \a + l) 0.24 + 1)3 U- j 
where g* is the number of relativistic degrees of freedom at the nucleation temperature. 
Assuming a detonation, the bubble wall velocity is given approximately by B38l 



. 1/V3+ Va 2 + 2a/3 
v b {a) = — . (8) 



In addition [33 



k(a) 



1 + 0.715a 



4 3a 
0.715a + —J — 
27 V 2 



(9) 



The peak frequency is 



f^ 52 * w iw){^){mV Hz - <10) 

Analysis of two-bubble collisions suggests that the spectrum rises as / 2 8 and / _1 8 , below 
and above the peak frequency, respectively [35]. 

The subject of the GW production from first order phase transitions is still a field of 
active research and of continuous developments. The authors of [39 J developed a different 
modeling of the problem. Instead of performing numerical simulations of colliding 
bubbles, they considered the bubble wall velocity as a random variable. Although in this 
approach the collisions are not formulated in a deterministic way, the advantage is that 
the envelope approximation in this case is not implemented. The spectrum is 



where Vf = {v\, - 1/ V3)/(l - 0&/ V5), and s = l/{vi, y3). The peak frequency is 



/»- mxir SIooG5vlraoJ ^ Hz ' <12) 

Away from the peak frequency, the spectrum is multiplied by a factor 2.5/,. 3 / (1 +0.5/ r 2 +/ r 4 ' 8 ), 
where f r = //(0.87/ c ' oll ). Finally, new numerical simulations have been done recently with 
multicolliding bubbles BOB suggesting that the GW spectrum decreases like f~ l rather 
than / _1 8 . This calculation gives 

\ 2 / a \ 2 j /I0" xl/3 
jS j U + l/ 0A2 + v 2 h \ g 



with peak frequency 

65 v 1 rr 5 1 

,HJU00GeVAl00/ 1.8 - 0.1u b + v 



a =^x ^ mi^m" 6 t^^h, (1 4) 



b 
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The spectrum (according to this calculation) rises as f 3 for frequencies below the /" and 
falls off as f~ l for frequencies larger than f". 

As we have mentioned, the second source of GWs during a first order phase transition 
can be turbulence. When bubbles collide, the plasma is stirred up and develops the 
characteristics of a fully developed turbulence. This means that a cascade of eddies is 
created in the plasma. Large eddies (of the size of the system or the stirring source) 
are formed and after a few revolutions, they break down to smaller eddies until their 
size becomes equal to the damping scale. Fluid (nonrelativistic) turbulence has been 
found experimentally to agree with Kolmogorov's stochastic description. This description 
has been implemented in calculations of GW production due to turbulence ||34l l4~Tl - l4"3f , 
although the fluid in this case is relativistic. Although in all these calculations, the 
Kolmogorov spectrum is used in order to model turbulence, there is a sort of different 
philosophy between for example |4"T|, and [|42|. In the former, GWs inherit directly the 



frequency of the eddies, while in the latter GWs inherit the wave number of the eddies. It is 
easy to see that the two approaches are not equivalent. The characteristic frequency of the 
eddies is co\ = v s /l (with v s being the velocity of the fluid in the eddy and I the characteristic 
length of it). In the first approach, the GWs peak at the frequency of the largest eddy 
which is co l = v s /L (with L being the size of the stirring source). On the other hand, 
the wave number of the stirring source is k ~ 1 /L. If the GWs inherit the wave number 
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instead of the frequency of the eddies (as in the second approach), the GW spectrum 
(because its dispersion relation is co = \k\) peaks at frequencies co = l/L £ v s /L = co L (since 
v s + 1). Here we follow the approach presented in [44| for the GW production due to first 
order phase transitions in the early universe. In this framework the GWs inherit the wave 
number spectrum. However, we have also checked that very similar results are obtained 
when adopting the framework discussed in Il43ll45l . The GW density for v s < 0.5 is 

f<2» J, 
2v s f p <f<f P , 



n turh h 2 = 6.7 x 1(T 6 A? 



'H»\ 2 /100\ 1/3 



i f r 8/3 

1(1) - A</ 



and for v, > 0.5 is 



Du.rhh 2 = 6.7 x 10- 6 vy b I -4- 



H»\ 2 /100\ 1/3 / 1 



P \8- 



Av 2 



where 



ka 



2 

(|) , fp<f<8vU P , 
4v 2 ^} , 8v*f p <f 



4/3 + ka 



The peak frequency f p is ||44 



^^ 8xl0 " 6 ^(l)(ira)(iir Hz - 



(15) 



(16) 



III. GRAVITATIONAL WAVES FROM MINIMAL WALKING TECHNICOLOR 



The new dynamical sector we consider, which underlies the Higgs mechanism, is an 
SU(2) technicolor gauge theory with two adjoint technifermions [2J. The two adjoint 
fermions may be written as 



Ql 



u a 



U a R , D a R , a = 1,2,3, 



(17) 



with a being the adjoint color index of SU(2). The left-handed fields are arranged in 
three doublets of the SU(2) L weak interactions in the standard fashion. The condensate 



is (UU + DD) which correctly breaks the electroweak symmetry. The model as described 
so far suffers from the Witten topological anomaly ||46| . However, this can easily be 
addressed by adding a new weakly charged fermionic doublet which is a technicolor 
singlet 0. 

In UU we constructed the effective theory for MWT including composite scalars and 
vector bosons, their self-interactions, and their interactions with the electroweak gauge 
fields and the SM fermions. We have also used the Weinberg modified sum rules [|47| to 
constrain the low energy effective theory This extension of the SM was thereby shown 
to pass the electroweak precision tests. Near the finite temperature phase transition the 
relevant degrees of freedom are the scalars and hence we will not consider the vector 
spectrum or that of the composite fermions. 

The relevant effective theory for the Higgs sector at the electroweak scale consists, in 
our model, of a composite Higgs and its pseudoscalar partner, as well as nine pseudoscalar 
Goldstone bosons and their scalar partners. These can be assembled in the matrix 

o + i© 



M 



+ V2(z'n fl + n fl ) x a 



E , (18) 



2 

which transforms under the full SLT(4) group according to 

M^uMu 7 , with ueSU(4). (19) 

The X a 's, a = 1, . . . , 9, are the generators of the SLT(4) group which do not leave the vacuum 
expectation value (VEV) of M invariant. 



The electroweak subgroup can be embedded in SU(4), as explained in detail in [48 1 
The new Higgs Lagrangian is 



Xniggs = ^Tr [D,M& l M f ] - <V{M) + £etc , (20) 
where the potential reads 

<V(M) = -^-Tr[MM + ] + ^Tr[MM + ] 2 + A'Tr[MM + MM + ] 

- 2A" [Det(M) + Det(M + )] , (21) 

and Xetc contains all terms which are generated by the ETC interactions, and not by the 
chiral symmetry breaking sector. 



10 



We explicitly break the SU(4) symmetry in order to provide mass to the Goldstone 
bosons which are not eaten by the weak gauge bosons. Assuming parity invariance, 

2 

Tit 

Xetc = -jp Tr [MBM f B + MM 1 ] + ■ • • , (22) 



where the ellipses represent possible higher dimensional operators, and B is a constant 
matrix that commutes with the SU(2) L xSU(2) R xU(l) v generators. 

The potential *V(M) is SU(4) invariant. It produces a VEV which parametrizes the 
techniquark condensate, and spontaneously breaks SU(4) to SO(4). In terms of the model 
parameters the VEV is 

" 2 = <"> 2 = ITa^' <23) 

while the Higgs mass is 

M 2 H = 2m 2 . (24) 

The linear combination A + A' - A" corresponds to the Higgs self-coupling in the SM. The 
three pseudoscalar mesons IT*, IT correspond to the three massless Goldstone bosons 
which are absorbed by the longitudinal degrees of freedom of the W ± and Z boson. The 
remaining six uneaten Goldstone bosons are technibaryons, and all acquire tree-level 
degenerate masses through (not yet specified) ETC interactions [59 j: 

K uu = M 2 UuD = M 2 Udd = ml TC . (25) 

The remaining scalar and pseudoscalar masses are 

M@ = 4v 2 A" 

M 2 A± =M 2 A0 = 2v 2 (A' + A") (26) 

for the technimesons, and 

Ml =M 2 - =Ml =m 2 ETC + 2v 2 (A' + A") , (27) 

fiuu n UD n DD E1> - 

for the technibaryons. Ref . ||49l provides further insight into some of these mass relations. 
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A. Effective potential for MWT 



The electroweak phase transition is studied by using the effective potential method. 
We include temperature dependent corrections of the effective potential up to one-loop 
level and ring resummation following Arnold & Espinosa [|5Ql|. We follow otherwise 
the conventions of [23] but use a slightly different method for estimating the temperature 
dependent one-loop correction, which involves a combination of the high temperature and 
low temperature asymptotic series ||5T1 (see Appendix [A] for details). This improves the 
potential at low temperatures. This is necessary here since strong GW production requires 
strong first order phase transitions, which typically means that the critical temperature 
is much smaller than the electroweak scale. Moreover, we use the actual nucleation 
temperature T», which can be considerably smaller than the critical temperature T c (where 
the symmetric phase and broken phase vacua are exactly degenerate) for strong first order 
transitions. 

We include in the analysis the heaviest standard model particles, the top quark and 
the weak gauge bosons. In addition, we consider the fourth family leptons, and a few 
composite scalar states that are made of techniquarks. The scalar states are expected to 
be the lightest states of the technicolor theory and have masses near the electroweak scale 
with strong coupling to the chiral condensate, which is identified with the expectation 
value of the composite Higgs. Hence they are the most prominent states for the dynamics 
of the electroweak phase transition. Of the two scenarios presented in [23J - light and 
heavy ETC masses - we only consider the latter one since it was seen to produce a stronger 
phase transition, potentially leading to stronger GWs. In this scenario the baryonic 
Goldstone bosons of the SLT(4) — > SO(4) chiral symmetry breaking are decoupled from 
the phase transition because of an ETC mass contribution that is much larger than the 
electroweak scale. The remaining eight scalar states include the composite Higgs a and 
its pseudoscalar partner as well as the Goldstone bosons n that are eaten by the gauge 
bosons, and their scalar partners A. 



The effective potential is then calculated as outlined in Appendix A and in |23"|]. We 
include the states listed above, except for the zero-temperature one-loop correction, where 
we leave out the (negligible but) infrared divergent contribution from the massless Gold- 
stone bosons. 
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FIG. 1: The nucleation temperature T* (left column) and the strength of the transition (p*/T* (right 
column) for MWT in the M H -M & plane for M Ar M f = 150 GeV and 350 GeV, as indicated in the 
labels. In the white regions the phase transition is either second order, very weakly first order, or 
does not occur at all. 
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B. Results 



Let us first comment on the parameter space. One of the couplings m 1 , A, A', and A" of 
the Higgs sector in the effective theory is fixed by the requirement v = 246 GeV, while the 
rest can be expressed in terms of the Higgs mass Mh, the mass M® of the particle, and 
the mass Ma of the scalar partners of the Goldstone bosons. Additional free parameters 
in our setup are the masses of the fourth family leptons which we assume to be equal and 
denote by M/. 

We follow [23J and present the results on the M H - M® plane while keeping M A 
and Mf fixed at reference values 150 GeV and 350 GeV. Fig. [I] shows our results for 
the nucleation temperature T» (left) and the strength of the phase transition (p*/T*. We 
compared these results to our earlier estimate [23j for the strength (p c /T c of the phase 
transition at the critical temperature, where the two vacua are exactly degenerate. There 
is a notable difference only when the first order transition is very strong, (£»/T* > 1. In this 
case typically <p*/T* > (j) c /T c . The region where the transition is strong enough to drive 
electroweak baryogenesis (<p/T > 1) is practically unchanged. Fig.[2]shows the results for 
the parameter a characterizing the produced latent heat (left), and the parameter j8/fi* 
characterizing the rapidness of the transition (right), which are specifically important 
for the production of GWs. It is seen that the values of the parameters are strongly 
correlated. As in other models (see, for example, I1281P , strong first order transition, with 
sizeable <p*/T*, generally means low /3/H* and large a, which are required for eminent 
production of GWs. The phase transition is at its strongest near a critical line on the 
Mh~M® plane, where (within our approach) T* — > and <p,/T, — > oo. However, values of 
a > 0.5 (and jS/H* - 50 . . . 1000 ) that are required [32] for the waves to be detectable at 
LISA B52| are obtained only in a very narrow slice of the parameter space. 

We observe that the value assumed by /S/H* increases substantially as the strength of 
the phase transition decreases. We find, for strong phase transitions, agreement with the 
general prediction for this ratio put forward in ||37l |53| while for weak phase transitions 
we find agreement with the thin wall approximation. 

In Fig. [3] we have plotted the gravitational spectrum of a MWT theory with a = 0.2, 
j6 /H» = 300 and T* = 60 GeV. This is an example of a strong first order phase transition 
with (p*/T* - 4, which can be obtained without fine-tuning the parameter values to be 
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FIG. 2: The parameters a (left column) and jS/H» (right column), which characterize the GW 
production, for MWT in the M H -M @ plane for M Ar M f = 150 GeV and 350 GeV, as indicated in the 
labels. 
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FIG. 3: The density of produced GWs pew = ^cwh 2 in the case of the MWT as a function of the 
frequency (in Hz). Dashed lines represent the expected sensitivity of LISA and BBO, while solid 
lines represent the gravitational spectrum of bubble collision and turbulence combined. The three 
solid lines from thinner to thicker correspond to the gravitational spectrum with bubble collisions 
given, respectively, by Eqs. j7|), |13) , ( [TT| ). The values of the parameters are given in the text. 

unnaturally close to the critical line. This set of numbers can be derived from the effective 
potential if we choose, for example, M A = M f = 150 GeV, M H 150 GeV, and M - 583 
GeV. For these values of M H , M A , and Mf the phase transition disappears for M & > 592 
GeV, so the chosen value of M@ lies within 10 GeV from the critical line. 

IV. GRAVITATIONAL WAVES FROM MULTIPLE PHASE TRANSITIONS 

An interesting scenario which may arise in strongly interacting extensions of the stan- 
dard model is that there can be several phase transitions at temperatures close to the 
electroweak scale ||24|. Only one of these transitions needs to break the electroweak sym- 
metry, while all of them can produce GWs, which can lead to a complex GW spectrum 
with several peaks from the various transitions. We shall consider here the case of Ultra 
Minimal Technicolor (UMT) iTTOll , where chiral symmetry breaking can proceed in two (or 
possibly three [24 j) steps, related to the two sectors of matter in this theory. 

In general, the two sectors of UMT can talk to each other, since the matter in the 
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different sectors interact via the strong technicolor dynamics. This interplay is realized 
at the effective Lagrangian level by terms which involve composite bound states from 
both sectors. To simplify the discussion on the production of GWs, we shall here omit 
these terms (set 6 = = 6' below) and assume that the two sectors are decoupled. Then 
the electroweak gauge symmetry feels only one of the sectors (the one with the order 
parameter 04 below), while the other (the one with the order parameter 02) is decoupled 
from the electroweak dynamics. In this case, the formalism for GW production presented 
above is directly applicable for UMT The effect of the interactions between the two sectors 
was studied in ||24ll25l| . The interactions were found to lead to a rich phase diagram with 
the possibility of breaking the electroweak symmetry twice (and restoring it once) as 
the universe cools down, while the phase transitions were typically weaker than in the 
decoupled case. Hence we expect that also the produced GWs are at their strongest in the 
scenario investigated here. However, generally even within UMT, there is the possibility 
of having three phase transitions and consequently a richer GW spectrum than the one 
presented here. 

Let us make a brief general comment on multiple (first order) transitions in the early 
universe. An important question is if the transitions can be treated separately or if 
there are bubbles related to different transitions present simultaneously. Recall from the 
discussion above that the nucleation probability density is ~ T 4 e~ S£ and consequently the 
time scale of nucleation is given by the inverse of /3 = —dS^/dt. The scale of change in 
temperature is thus 

AT T .tt (28) 

j6 dt j6 v ' 

Since we will have jS/H* > 100 (typically even /3/H* > 1000) and T* is around a few hundred 
GeV, AT is less than a few GeV. Hence the transitions can be practically treated as separate, 
unless the underlying dynamics forces them to be exactly simultaneous [24, 25J. Notice 
that the nucleation probability depends exponentially on Se', therefore the transition will 
end very quickly when the temperature difference with respect to the start of nucleation 
exceeds T»H*/j6. 



17 



A. Effective theory for UMT 



The model proposed in IITOl consists of an SU(2) gauge group with two Dirac fermions 
belonging to the fundamental representation and two Weyl fermions belonging to the 
adjoint representation. In order not to be in conflict with the Electroweak Precision Tests 
only the fundamental fermions are charged under the electroweak symmetry 

We shall only consider the effect of the composite scalar mesons which are expected to 
be the lightest particles in the theory Their masses have the strongest dependence on the 
vacuum expectations values of the Higgs fields. 

The relevant degrees of freedom are efficiently collected in two distinct matrices, M4 
and M2, which transform as M4 — > g 4 M 4 gJ and M 2 — > giMig T 1 with g 4 £ StT(4) and 
g2 e SU(2). Both M 4 and M 2 consist of a composite isoscalar and its pseudoscalar partner 
together with the Goldstone bosons and their scalar partners 



The notation is such that X4 and X 2 are the broken generators of StT(4) and SU(2) 
respectively. Also cr 4 and ©4 are the composite Higgs and its pseudoscalar partner while 



simply substitutes the index 4 with the index 2. 

To describe the interaction with the weak gauge bosons we embed the electroweak 
gauge group in StT(4) as done in ||48|. Because of the choice of the electroweak embedding 
the weak interactions explicitly reduce the SLT(4) symmetry to SU(2) L x U(l)y x U(1)tb 
which is further broken to li(l) em X U(1) T b via the technicolor interactions. U(1) T b is the 
technibaryon number related to the fundamental fermions. The remaining SU(2) X IT(1) 
spontaneously breaks, via the extra technifermion condensates, to SO(2) x Z 2 . Here 
SO(2) = U(l) is the technibaryon number related to the adjoint fermions. 

We are now in a position to write down the effective Lagrangian. It contains the kinetic 
terms and a potential term 





i = 1,2. 



(29) 



(30) 



and are the Goldstone bosons and their associated scalar partners. For SU(2) one 



X = -Tr[D^M 4 D^] + -Tr[^M 2 <9^]-T(M 4 ,M 2 ) , 



(31) 
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where the potential is 



<V (M 4 , M 2 ) = - ^Tr [M 4 Ml] + ^Tr [M 4 M+] Z + A^Tr [M 4 M+M 4 Ml] 



- y Tr [m 2 MI] + ^Tr [m 2 M+] Z + A 2 Tr [a^M+A^M+J (32) 

+ ^Tr [M 4 M+] Tr [m 2 M+] + 46' [(detM 2 ) 2 Pf M 4 + h.c] . 

We shall from now on set 5 = = 5' so that the two sectors are decoupled. 

Once M 4 develops a vacuum expectation value the electroweak symmetry breaks and 
three of the eight Goldstone bosons - n°, n + , and II" - will be eaten by the massive gauge 
bosons. In terms of the parameters of the theory the vacuum states <cr 4 ) = c 4 and (a 2 ) = v 2 
which minimize the potential are 

m\ = (A. + A;)^, (33) 
m\ = (A 2 + 2A' 2 )v 2 2 . (34) 

For the model to be phenomenologically viable some of the Goldstone bosons must 
acquire a mass. Here we parametrize the ETC interactions by adding at the effective 
Lagrangian level the operators needed to give the unwanted Goldstone bosons an explicit 
mass term. 

The effective ETC Lagrangian breaks the global SLT(4) x SU(2) x U(l) symmetry. To 
construct the required ETC terms at the effective Lagrangian level, we find it useful to 
split M 4 (M 2 ) - form invariant under LT(4) (U(2)) - as follows: 



A 



M4-M4 + iPi , and M 2 = M 2 + zP 2 , 



with 



M 4 = 
M 2 = 



V2 



E 4/ P 4 = 
E 2/ P 2 = 







V2 



E 4 , 
E 2/ 



/ = 1,...,5, 
1 = 1,2. 



(35) 

(36) 
(37) 



M 4 (M 2 ) as well as P 4 (P 2 ) are separately SLT(4) (SLT(2)) form invariant. A set of operators 
able to give masses to the electroweak neutral Goldstone bosons is 



m 



4,ETC 



ETC - 



Tr [M 4 B 4 M+B 4 + M 4 M+] 



m 



2,ETC , 



-m\ ETC [Pf P 4 + Pf P+] 



m 



1,ETC 



Tr [m 2 B 2 M\B 2 + M 2 M+] 
[det(P 2 ) + det(P+)] , (38) 
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where B4 is the diagonal SLT(4) generator that commutes with the electroweak generators, 
and B 2 = t 3 is the diagonal generator of SU(2). The masses of the two Higgs particles 

M 2 Hi = 2m\ = 2(\, + A\)v\ (39) 
M 2 Hi = 2m\ = 2 (A 2 + 2A' 2 ) v\ (40) 

are unaffected by the addition of the ETC low energy operators. The rest of the spectrum 
is: 

M k D = m l,ETC > K M = m l,ETC > M k = KeTC = K 2 ( 41 ) 

for the pseudoscalar partners and the Goldstone bosons that are not eaten by the massive 
vector bosons and: 

M L = M no = M 2 n± = 2 (A' 4 vl + 5'vf) + m\ ETC , (42) 
M 2 „ aa = Av\ (A' 2 + b'v\) + m{ ETC , (43) 

for the scalar partners. 

B. Effective potential analysis for UMT 

For the effective potential in UMT we use the same approach as outlined above for 
MWT: we include one-loop corrections with ring resummation for the bosonic degrees of 



freedom as presented in Appendix A and in [25 J. The particle spectrum includes the top 
quark, the weak gauge bosons, and the lowest scalar states of the theory presented above. 
For UMT the ETC mass scale must be of the order of the electroweak scale in order to 
obtain first order transitions. 

The effective Lagrangian of scalar particles in UMT includes several free parameters. 
However, since we study only the case where the potentials related to the fundamental 
and adjoint techniquarks are decoupled, we have 6 = = 6' and each of the transitions 
depends only on a certain subset of parameters. The potential of the SLT(4) sector (funda- 
mental quarks) is characterized by the Higgs mass Mh 4 , the masses of the scalar partners 
of the Goldstone bosons, and the ETC masses. Since this sector breaks the electroweak 
symmetry, the zero-temperature VEV U4 must equal the electroweak scale 246 GeV. The 
dependence on the ETC masses is relatively weak, and we fix all of them to be 150 GeV. 
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FIG. 4: The nucleation temperature T», the strength of the phase transition <p*/T*, the parameter 
a, and the parameter jS/H» (from top to bottom) in the Mh-AMu plane for the "4 transition" (left) 
and for the "2 transition" (right) of UMT. We fixed ETC masses at 150 GeV and used v 2 = 300 GeV. 
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For the (dynamical contribution to the) masses of the scalar partners of the Goldstone 
bosons we use 

(AM n J 2 = M 2 nuo - m\ ETC = 2X\v\ . (44) 

Similar parameters characterize the potential of the SU(2) sector. In addition to the Higgs 
mass Mh 2 and 

(AMn 2 ) 2 = Mjj u - m\ ETC = W 2 vj (45) 

the zero-temperature value of the condensate is now a free parameter (while V4 was 
fixed to 246 GeV). Notice also that since only the fundamental techniquarks are charged 
under the electro weak gauge symmetry, the SU(2) — > SO(2) transition is independent of 
the standard model parameters. 

C. Results 

The results for UMT are shown in Fig. |4j The rows from top to bottom show the 
behavior of the nucleation temperature T», the ratio <p*/T* at the nucleation tempera- 
ture, the parameter a characterizing the produced latent heat, and the parameter j8/H* 
characterizing the rapidness of the transition, respectively. The left-hand plots give the 
parameters for the SLT(4) — > Sp(4) transition (the one coupled to the electroweak), while 
the right-hand plots are for the SU(2) — > SO(2) transition. 

In general the plots are very similar as in MWT above. The main difference is that 
both of the transitions here are weaker than in MWT: in particular, the crucial parameter 
a < 0.02. Therefore, it seems unlikely that detectable GWs could be produced in the UMT 
model. The maximal value of the a parameter can be slightly enhanced by optimizing 
the choice for the ETC masses, possibly by adding new ETC operators, and by increasing 
the value of V2 (which was fixed to 300 GeV above in Fig. |4|). 

Despite the weakness of the transition, there is however a very intriguing scenario. 
UMT admits successive phase transitions that can be of first order and occur at different 
temperatures. This in principle means that the gravitational spectra of the two (or more 
in general) phase transitions can peak at well separated frequencies. If the parameters a 
and jS are such that the spectrum of one phase transition does not completely cover the 
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FIG. 5: The density of produced GWs pew = ^cwh 2 in the case of the UMT as a function of 
the frequency (in Hz). Dashed lines represent the expected sensitivity of LISA and BBO, while 
solid lines represent the gravitational spectrum of bubble collision and turbulence combined. The 
thin and thick solid lines correspond to the gravitational spectrum with bubble collisions given, 



respectively by Eqs. ( |13[ l, and (111. The spectrum from the 04 (electro weak) transition peaks at 
roughly / = 0.05 Hz, while the spectrum from the 02 transition peaks around / = 0.005 Hz. The 
values of the parameters are given in the text. 



one of the second, then multiple peaks can be potentially seen. Such a case is depicted in 
Fig. |5j In this case, one phase transition, i.e., the one associated to 04 (the one breaking 
the electroweak symmetry), has a = 0.025, p7H* = 3000, and T» = 120 GeV, while the 
second phase transition associated to 02 has a = 0.005, jS/H* = 200, and T* = 200 GeV. 
This set of values can be deduced from the effective potential with M Hi 130 GeV and 
AMn 4 - 350 GeV for the "4 sector" and with M Hl « 180 GeV and AM Ui - 635 GeV for 
the "2 sector." The ETC masses were 150 GeV and = 300 GeV as above. As seen from 
Fig. [4] this is a rather optimal scenario with a low Higgs mass Mh 4 and parameters for 
the 02 transition rather near the critical line. In Fig. [5] we see the existence of multiple 
peaks that span almost two orders of magnitude in frequency (from 10" 3 to 10" 1 ). The 
peaks at higher frequency are due to the transition in the 04 sector, which also breaks 
the electroweak symmetry, while the lower peaks are produced by the 02 transition. 
Unfortunately as we see, the spectrum lies below the expected sensitivity of BBO. As 
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mentioned above, the spectrum can be enhanced by choosing optimal values for the ETC 
masses and the vacuum expectation value v 2 . However, in all the configurations which 
we have checked, this does not change the conclusion: the nontrivial structure of the 
spectrum is hardly visible at BBO. Nevertheless, this is an interesting case for two reasons 
(apart from hoping for a better sensitivity of BBO). First, one can consider an underlying 
theory which has a larger number flavors than UMT. The presence of extra fermions in the 
theory can strengthen the phase transitions. In such a case, multiple peaks can be above 
the BBO sensitivity. Second, the estimation of GWs from first order phase transitions 
is far from conclusive. For example, apart from bubble collision and pure turbulence 
as we considered here, there is a possibility of producing GWs via primordial magnetic 
fields ||54] - l56| , with potentially larger amplitudes. 

V. CONCLUSION 

In this paper we studied the production of GWs from first order phase transitions of 
theories that dynamically break the electroweak symmetry. Although our setup is general, 
we looked, in particular, at two different models, i.e., the Minimal Walking Technicolor, 
and the Ultra Minimal Technicolor, that have been studied extensively. For MWT we 
found that there is parameter space, at the effective Lagrangian level, for a sufficiently 
strong first order phase transition that produces GWs detectable at BBO. In this case, 
however, the phase transition is so strong that if a baryogenesis mechanism takes place, 
sphalerons would not be able to wash out the produced asymmetry. First principle lattice 
simulations will be able to disentangle, in the near future, the order of the phase transition 
of the underlying gauge theory. We should also stress that the low energy effective theory 
used here for the MWT Lagrangian can also describe the low energy effective theory for an 
SO(4) gauge theory with 2 Dirac fermions in the fundamental representation of the gauge 
group. The latter gauge theory has two advantages over the traditional SU(2) gauge 
theory with fermions in the adjoint representation. It is not expected to be conformal 
||57|| and does not feature technigluons-techniquark bound states ||22| with potentially 
dangerous fractionally charged states. We have also discovered that, for reasonable 
values of the parameters, UMT seems not to be able to provide a very strong first order 
phase transition. However, UMT undergoes successive phase transitions, which can 
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produce a gravitational spectrum of multiple peaks spanning two orders of magnitude 
in the frequency. This provides a very characteristic signal that can differentiate strongly 
coupled theories with multiple first order phase transitions. 
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Appendix A: The effective Potential 



The effective potential is obtained by adding to the tree-level potential O) the one-loop 
correction 



V(o) = V^(o) + V%(o) + V?(o, T) . 



(Al) 



For brevity, we denote by a the expectation values of the Higgs field(s) that characterize 
the techniquark condensate(s). For MWT we have a single condensate o, while for UMT 
we identify a = {04, 02} where cr, refer to the condensates of the two sectors in UMT The 
standard zero-temperature one-loop contribution to the potential reads 



(A2) 



where the index i runs over all of the mass eigenstates and n, is the multiplicity factor for 
a given scalar particle n\,, while for Dirac fermions it is -4 times the multiplicity factor of 
the specific fermion tif. The function fi is 



log 



M?(g) 3 
Mf(v) 2 



+ 2M 2 (a)M 2 (v), 



(A3) 



where M 2 (o) is the background dependent mass term of the z'th particle and ot=o = v. 

The one-loop, ring-improved, correction can be divided into fermionic, scalar, and 
vector contributions, 



V? = v? f + V?> + V? 

T T f T b T 



gauge 



(A4) 



We use the extrapolation method introduced in [pTH for evaluating the one-loop correction. 
At high temperatures we expand in M//T, which gives for the fermions 



16n 2 



f 



log 



M){o) 
j, 2 



c f 



? V » M 2 (nYT 2 V (~ M f (a) ) gj - WW ~ V M-i A , A « 
-2 ^ n f M f (a)T ^ (7mm + v ( 2 ~ l ) > < A5 > 



(2/)!!(Z + l) 



/ '=2 

where ct — 2.63505. For the bosons (including the electroweak gauge bosons, for which 
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fib = 3) we write instead 

b 

b 

M^T 2 V- /-^(^) V (21 - 3)!!C(2Z - 1) 



log - c b 



Hb 2 lu\4n^ (2/)!!(Z + l) ' (A5) 



+ 

b ~ 1=2 



where C\, — 5.40762. At low temperatures we use for both bosons and fermions the 
asymptotic expansion of the one-loop correction 

The extrapolated one-loop correction reads for each fermion 

V*\ = &L- J V^\ h (N = 4) + 40 - Xf \ (V^\(N = 3)-5 f ), (A8) 

where is the step function. The parameter Xf — 2.21605 and the small correction 5f - 
-7.90454 x 10" 4 were fixed by requiring the function to be continuous and differentiable 
with respect to Mf at M^/T 2 = Xf. For the scalars we resum the contribution of the ring 
diagrams. Following Arnold and Espinosa ||50l we write 

V» = O (*> - V» bt) (N = 3) + - *) ( V»CN = 3) - 6„) 

+^(M>)-M>,T)) , (A9) 

where x\, - 9.47134, b\, 3.1931 x 10" 4 , and Mb(o, T) is the thermal mass which follows 
from the tree-level plus one-loop thermal contribution to the potential. For the gauge 
bosons we set n b = 3: 

vf %b - e(, t - W)v<V(N = 3 ) + e(W-, t )(v?>(N = 3)-a t ) 

+ is( M >'- M U' J ' T >) ■ (A10 > 

Here Mi ig b(o, T) is the longitudinal mass of a given gauge boson and we have Mi ig b(o, T = 
0) = M g b(o) while the transverse mass receives only a suppressed temperature dependent 
correction which we have neglected. 
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For a more complete presentation with explicit expressions for the masses see Il23ll25l . 
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